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I. INTRODUCTION 



The surprising discovery of new states decaying to 
J/yjir + n~ [1, 2] has renewed interest in the field of 
charmonium spectroscopy, since the new resonances are 
not easy to accommodate in the quark model. In par- 
ticular, the BaBar experiment discovered a new broad 
state, F(4260), decaying to J/ipTr + TT~ in the initial-state- 
radiation (ISR) reaction e+e" -> i ISR Y (4260). The 
quantum numbers J PC = 1~~ are inferred from the sin- 
gle virtual-photon production mechanism. Further struc- 
tures at 4.36 GeV/c 2 [3, 4] and 4.66 GeV/c 2 [4] have 
been observed in the yj(2S)ir + n~ mass distribution from 
the reaction e + e~ — ► "fjsR'4'(2S)'K + Tr^. Charmonium 
states at these masses would be expected [5, 6] to decay 
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predominantly to DD, D*D, or D*D* [7]. It is pecu- 
liar that the decay rate to the hidden charm final state 
J/V>7r + 7r~ is much larger for the V(4260) than for excited 
charmonium states [8], and that at the V(4260) mass 
the cross section for e + e~ — > hadrons exhibits a local 
minimum [9]. Several theoretical interpretations for the 
y(4260) have been proposed, including unconventional 
scenarios: quark-antiquark gluon hybrids [10], baryo- 
nium [11], tetraquarks [12], and hadronic molecules [13]. 
For a discussion and a list of references see, for example, 
Ref. [14]. 

This work explores ISR production of DD, D*D, and 
D*D* final states for evidence of charmonium states and 
unconventional structures. This follows an earlier BABAR 
measurement of the DD cross section [15]. A study by 
the Belle collaboration of the DD, D*D, and D*D* final 
states can be found in Rcfs. [16, 17]. Recent measure- 
ments of the e + e~ cross sections can be found in Ref. [18]. 

We also measure for the first time branching fractions 
of high mass charmonium states, other than V(4260), 
for which little information exists [9], and compare our 
measurements with theoretical expectations [5, 6, 14]. 

This paper is organized as follows. In Section II we 
give a short description of the BABAR experiment and in 
Section III we describe the data selection. Section IV 
is devoted to the selection of the D*D final state and in 
Section V, we present the mass resolution, reconstruction 
efficiency, and measured cross sections. In Section VI we 
describe the D*D* cross section measurement while in 
Section VII we present the DD data. The description of 
the fit of the three channels is described in Section VIII, 
while Section IX is devoted to the measurements of the 
ratios of branching fractions. Finally, in Section X, we 
compute the limit on production of Y"(4260) decaying to 
D*D and D*D* , and summarize conclusions in Section 
XL 



II. THE BABAR EXPERIMENT 

This analysis is based on a 384 fb _1 data sample 
recorded at the T(45) resonance and 40 MeV below 
the resonance by the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage rings. The BABAR de- 
tector is described in detail elsewhere [19]. We mention 
here only the parts of the detector which are used in 
the present analysis. Charged particles are detected and 
their momenta measured with a combination of a cylin- 
drical drift chamber (DCH) and a silicon vertex tracker 
(SVT), both operating within a 1.5 T magnetic field of a 
superconducting solenoid. The information from a ring- 
imaging Cherenkov detector combined with energy-loss 
measurements in the SVT and DCH provide identifi- 
cation of charged kaon and pion candidates. Photon 
energies are measured with a CsI(Tl) electromagnetic 
calorimeter. 



III. DATA SELECTION 

DD candidates are reconstructed in the seven final 
states listed in Table I. The D*° -> D°n° and £>*°_-> 
L>°7 decay modes are used to form D*°D° and D*°D*° 
candidates. The D*+ -» D°tt+ and D*+ -> D+n° de- 
cay modes are used to form D* + D~ and D* + D*~ can- 
didates. Table II summarizes the full decay chains used 
to reconstruct the D*D and D*D* candidates. 

For all final states, events are retained if the number of 
well-measured charged tracks, having a minimum trans- 
verse momentum of 0.1 GcV/c, is exactly equal to the 
total number of charged daughter particles. Photons are 
identified as EMC clusters that do not have a spatial 
match with a charged track, and that have a minimum 
energy of 30 MeV. Neutral pion candidates are formed 
from pairs of photons kinematically fitted with the n° 
mass constraint. K® candidates are reconstructed, with 
a vertex fit, in the tt + tt~ decay mode. The tracks corre- 
sponding to the charged daughters of each D candidate 
are constrained to come from a common vertex. Addi- 
tionally, for the D a -> K~ ~tt + tt channel, the D mass 
constraint is included in the fit, and for the D~ — > K®ir~ 
channel, a K® mass constraint is imposed. Reconstructed 
D candidates with a \ 2 fit probability greater than 0.1% 
are retained. Each DD pair is refit to a common ver- 
tex with the constraint that the pair originates from the 
e+e~ interaction region. Only candidates with a x 2 fit 
probability greater than 0.1% are retained. Background 
7r° candidates from random combinations of photons and 
other background channels are suppressed by requiring 
no more than one tt° candidate other than those at- 
tributed to the D° and D* decays. Similarly, we require 
in the event no more than one extra photon candidate, 
having a minimum energy of 100 MeV, apart from any 
photon attributed to D* or it decays. 

For D decay modes without a ir° daughter, the D- 
candidate momentum is determined from the summed 
three-momenta of the decay particles and its energy is 
computed using the nominal D mass value [9]. For 
the D° —> K~iT + ir channel, the 4-momcntum from the 
mass-constrained fit is used. Similarly, the D* momen- 
tum is determined from the summed three-momenta of 
the decay particles and its energy is computed using the 
nominal D* mass. 

The ISR photon is preferentially emitted at small an- 
gles with respect to the beam axis, and escapes detec- 
tion in the majority of ISR events. Consequently, the 
ISR photon is treated as a missing particle. We define 
the squared mass (M r 2 cc ) recoiling against the DD, D*D, 
and D*D* systems using the four-momenta of the beam 
particles (p e ±) and of the reconstructed D (po) and D* 

M rec = (Pe- +Pe+ -PDM ~ Pet')) 2 (1) 

This quantity should peak near zero for both ISR events 
and for exclusive production of e + e~ — > D*D or e + e~ — > 
D*D*. In exclusive production the D*D and D*D* mass 
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TABLE I: List of the reconstructed DD final states. 



N Channel First D decay mode Second D decay mode 



1 


D°D° 


D o _> K - n + 


27° 


-> K+%- 


2 


D°D° 




D° 


— ► tf+Tr"^ 


3 


D°D° 




D° 


— ► A + 7T~7r + 7r~ 


4 


D°D° 


Z5° - K-ir+n 


D° 


— ► K + -K~ir + -iv~ 


5 


D + D- 


D + -f 1^71-+ tt+ 


D~ 


— > A + 7r~7r~ 


6 


D+D- 




D~ 


-> K+X-tt- 


7 


D+D- 




D~ 





TABLE II: List of the D* D and D* D* reconstructed final states. The reconstructed D° decay modes are listed in Table I for 
the D*° D and D* D*° states. The column headed "Veto" lists ambiguities with the indicated channels, "Removed" indicates 
the fraction of events removed by the veto. 



N Channel 


First decay mode 






Second decay 


mode 




Veto 


Removed % 


8 D*°D° 


7T° -> D°7 












9-12 


5.9 


9 D*°D° 


D*° -> D ( V> 












11,12 


3.2 


10 D*°D*° 


Z>*° -> L>° 7 






D*° _D° 7 






9,11 


1.1 


11 D*°D*° 


D*° -> D%° 






D*° - D° 7 






8,10 


0.7 


12 D*°D*° 


Z>*° -> D%° 
















13 D*+D~ 


D*+ ^ flO T +, DO 


-f AT" 




£)- -> A-+7T- 










14 D*+D- 


D*+ -> D+tt , D+ 


-»■ K 


_ 7T + 7r + 


£)- -> A-+7T- 


7T~ 








15 D*+D~ 


D*+ -^D°tt+, D° 


-f A" 


A+ 


£T -> A+tt" 


7T~ 








16 D*+D- 




-» A" 


^+ 


£)- -> A+A 










17 D*+L>- 


£>*+ ^ D+7T , £,+ 


-f A" 


~TT + n + 


D~ -> K + K~ 










18 D*+D*~ 




-»■ if" 


^+ 


D*~ -> 27V 




A+TT" 






19 D* + D*~ 


75*+ ^ D+ 7r Q ) ^+ 


-> A" 


-7T+7T + 






A+TT" 







distributions peak at the kinematic limit. Therefore 
we select ISR candidates by requiring DD, D*D and 
D*D* invariant masses below 6 GeV/c 2 and |M r 2 ec | < 
1 GcV 2 /c 4 . 

We select D and D* candidates based on the candidate 
D mass and the mass difference Am = Mrj>» — Mb- The 
D and D* parameters are obtained by fitting the relevant 
mass spectra (see Fig. 1 for some Am distributions) using 
a polynomial for the background and a single Gaussian 
for the signal. Events are selected within ±2.5cx from the 
fitted central values, where a is the Gaussian width. For 
D* + — ► D°tt + , the selection criterion has been extended 
to ±6cr due to the presence of non-Gaussian tails. 

Because of our tolerance of an extra 7r° and/or 7, an 
ambiguity can occur for channels involving a D*° which 
is handled as follows. Each combination is considered as 
a possible candidate for channels 8-12 and D°D°. Monte 
Carlo simulations weighted by the DD, D*D, and D*D* 
measured cross sections [15-17] and branching fractions 
are used to optimize the selection criteria and estimate 
the fecdthrough of one channel to the other. A candidate 
is rejected if (a) it satisfies all the selection criteria for 
an ambiguous channel and (b) this rejection does not 
produce any significant loss in the channel under study 
and therefore can be classified as background. The list of 
channels rejected in case of ambiguities are indicated in 



the "Veto" column in Table II. The table also lists the 
fraction of events removed by these cuts in the |M 2 oc < 
1 GcV 2 /c 4 region. 

In the case of multiple D*° candidates, such as D*°D*° 
with both D*° — ► D°j, the candidate with m(D°j) clos- 
est to the nominal D*° mass is accepted. The charged 
D*D and D*D* modes, also listed in Table II, do not 
require such a procedure because backgrounds are negli- 
gible. 

IV. STUDY OF THE D* D FINAL STATE 

Figure 2 shows the D*D M 2 oc distributions after all 
the cuts for (a) D*°D°,D*° -> L>° 7 , (b) D*°D°,D*° -> 
D°7r°, and (c) D* + D~ . Clear peaks centered at zero 
with little background are observed, providing evidence 
of an ISR process. The number of background events 
in the |M 2 ec | < 1 GeV 2 /c 4 is estimated by fitting the 
M 2 ec distribution for each channel. The fits are per- 
formed using a 2 nd order polynomial for the background 
and a signal M r 2 c lineshape obtained from Monte Carlo 
simulations corresponding to the relative composition of 
the data. The event yields are obtained by subtract- 
ing the fitted backgrounds and integrating the result- 
ing M 2 CC distributions in the |M 2 CC | < 1 GeV 2 /c 4 re- 
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0.2 

Am(DV) (GeV/c 2 ) 



0.12 0.14 0.16 

Am(DV) (GeV/c 2 ) 




0.14 0.15 

Am(DV) (GeV/c 2 



0.16 0.12 0.14 0.16 

Am(DV) (GeV/c 2 ) 



TABLE III: Number of ISR candidates and purities for the 
different channels calculated in the range |M r 2 cc | < 1 GeV 2 /c 4 . 
The last column gives the value of the average efficiency e s 
at a mass of 4.5 GeV/c 2 . 



Channel 


bignal+ 
Background 


rurity(yoJ 


e x 1U 


D°D° 
D+D~ 


654 
199 


74.3 ± 1.7 

88.4 ± 2.3 




Total DD 


853 


77.6 ± 1.4 


25 


D*°D°, D*° -> D°-y 
D*°D° D*° -> D°n° 


460 
422 


75.4 ± 2.0 
84.4 ± 1.8 




Total D*°D° 


882 


79.7 ± 1.4 


4 


D*+D~ 


228 


100±° 


5 


Total D'D 


1110 


83.6 ± 1.1 




D*°D*° 


293 


69.3 ± 3.7 




D *+ D *- 


33 


100±° 




Total D*D* 


326 


72.1 ± 2.5 


1 



MASS RESOLUTION, EFFICIENCY, AND 
D* D CROSS SECTION 



FIG. 1: Am distributions for D* D candidates after applying 
the |M r 2 cc | < 1 GeV 2 /c 4 and m(D*D) < 6 GeV/c 2 selections, 
for (a) D*° -> D°j, (b) D*° -> D°ir°, (c) D*+ D°7r+ with 
D° -> if"7r + , and (d) D*+ -» D+tt with D+ -> A"-7r + 7r+. 
The shaded regions indicate the ranges used to select the D* 
candidates. 



gion. The resulting yields and fitted purities, defined 
as P = A^ignai/ (Signal + ^background), for each channel 
are summarized in Table III. 

The purity of each reconstructed D* channel is also 
demonstrated in Fig. 1, where the Am distribution is 
shown for D*D candidates with |M 2 CC < 1 GeV 2 /c 4 and 
D*D masses below 6 GeV/c 2 . The final selection of the 
ISR candidates is performed applying the Am selection 
criteria described above. 

The D*°D° mass spectrum is shown in Fig. 3(a), and 
the D* + D~ mass spectrum is shown in Fig. 3(b). Both 
spectra show an enhancement near threshold due to the 
presence of the ^(4040) resonance. 

The background shape for D*°D° candidates is ex- 
plored using the M r 2 cc sideband region, 1.5 < M 2 ec < 
3.5 GeV 2 /c . The D*°D° mass spectrum for these events, 
normalized to the background estimated from the fit to 
the M r 2 ec distribution, is presented as the shaded his- 
togram in Fig. 3(a). This background has been fitted 
with a threshold function: 



id f ~\ /™ ~~ \a-\-0m —ym — Sm 2 — em 3 

B[m) = [m — m t h) e 



(2) 



where m th is the threshold D*°D° mass. The D*+D~ 
final state is consistent with having zero background. 



In order to measure efficiencies and D*D mass resolu- 
tions, ISR events are simulated at five different values of 
the D*D invariant masses between 4.25 and 6.25 GeV/c 2 . 
These events are simulated using the GEANT4 detector 
simulation package [20] and are processed through the 
same reconstruction and analysis chain as real events. 
The mass resolution is determined from the difference 
between generated and reconstructed D*D masses. The 
D*D mass resolutions are similar for all channels and 
increase with D*D mass from 5 to 10 MeV/c 2 . 

The mass-dependent reconstruction efficiency for chan- 
nel i, ei(m D «jj), is parameterized by a second-order poly- 
nomial, and is multiplied by the decay branching fraction 
B, [9] to define 



(3) 



These values are weighted by Ni(m D , the number of 
D*D candidates in channel i, to compute the average 
efficiency as a function of m D ,j=j, 



"D'D 



^=1 ef(m D , ff ) 



(4) 



where n is the number of decay modes. In this case 
we have eight D*°D° chanels (1-4 with D*° -> D°j 
and D*° -> L>°7T ) and two D*+D~ channels (13, 14). 
Representative values of e B , computed at a mass of 4.5 
GeV/c 2 , are displayed in Table III. The sample sizes for 
the Cabibbo-suppressed decay modes (15, 16, and 17 in 
Table II) are very small (32 events total) and comprise 
14% of the D* + D~ sample. The efficiency for these de- 
cay channels has not been directly computed; instead, 
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-5 5 -5 5 -5 5 

M 2 rec (GeV 2 /c 4 ) M 2 rec (GeV 2 /c 4 ) M 2 rec (GeV 2 /c 4 ) 



FIG. 2: Distribution of M r 2 oc , the mass recoiling against the D'D system, for (a) D*° — ► D°j, (b) D*° — ► D°ty°, and (c) 
D* + D~ candidates. The curves are the results from the fits described in the text. 



these modes arc treated as having the mean efficiency 
of the Cabibbo-allowed channels 13 and 14. The ten 
D'D channels, after correcting for efficiency and branch- 
ing fractions, have yields that are consistent within the 
statistical errors. 

The D*D cross section is computed using 

dN/dm Dt jj 



V e +e-->D*D( m D*D) = 



(5) 



e B {m D , D )dC/dm D , D 

where dN/dm D ,^ is the background-subtracted yield. 
The differential luminosity is computed as [21] 

dC 2m 



D* D 



dm 



D'D 



TTX 



■{Ms/ml) -l)(2-2x + x 2 ), (6) 



where s is the square of the e+e~ center-of-mass energy, 
a is the fine-structure constant, x = 1 — m 2 Dr -/s, m e is 
the electron mass, and £ is the integrated luminosity of 
384 fhT 1 . The_cross sections for D*°D°, D'+D', and 
combined D*°D° and D* + D~ are shown in Fig. 4. A 
clear V(4040) resonance is observed. 

The systematic uncertainties on the cross sections, 
10.9% for D*°D° and 9.3% for D*+D-, include un- 
certainties for particle identification, tracking, photon 
and 7T° reconstruction efficiencies, background estimates, 
branching fractions, and a potential inaccuracy in the 
simulation of extraneous tracks, photon and n° candi- 
dates. The uncertainty due to the ISR selection has 
been estimated by narrowing the M r 2 c allowed range to 
0.7 GeV 2 /c 4 . All contributions are added in quadrature. 
Systematic uncertainties are summarized in Table IV. 

The D*°D° and D* + D~ cross sections have similar 
features and consistent yields. Integrating the cross sec- 
tions from threshold to 6 GeV/c 2 , we obtain 



a(D*+D~) 
<7(D* a D°) 



0.95±0.09 stot ±0.10 SJ/st , (7) 



TABLE IV: Systematic errors, given as fractional errors ex- 
pressed in %, in the evaluation of the D'D cross section. 



Effect 


D*°D° 


D*+D~ 


Background subtraction 


2.6 


3.0 


Branching fractions 


7.4 


4.6 


M r 2 ec cut 


2.2 


0.0 


Particle identification 


1.8 


2.1 


Tracking efficiency 


2.2 


3.3 


Extraneous tracks 


5.7 


5.7 


7r and 7 reconstruction efficiency 


3.4 


3.0 


Extraneous tv° and 7 


0.5 


0.8 


Total 


10.9 


9.3 



consistent with unity. In this calculation systematic er- 
rors related to the M r 2 ec selection criteria and tracking 
efficiency have been ignored because they largely cancel 
in the ratio. 



VI. STUDY OF THE D'D* SYSTEM 

A similar analysis is carried out for D*D* channels. 
Figure 5 shows the Am distributions for D*D* candi- 
dates with |M r 2 cc | < 1 GeV 2 /c 4 and D*D* masses below 
6 GeV/c 2 . The peak at threshold in Fig. 5(a) is due to 
background from D*° — > D°ir° where one 7 from the low 
momentum ir° is lost. 

We select the two D* candidates and reject candidates 
reconstructed in any of the modes listed in the "veto" 
column in Table II. Figure 6 shows the D*°D*° M 2 oc 
distributions for channels 10-12. 

The total D*°D*° and D*+D*- M 2 oc distributions are 
shown in Fig. 7. The number of background events for 
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FIG. 3: (a) D*°D° and (b) D* + D~ mass spectra. The shaded 
histogram in (a) is obtained from the M r 2 cc sideband region 
1.5 < M r 2 cc < 3.5 GeV 2 /c 4 normalized to the background 
estimated from the fit to the M r 2 oc distribution. The curve is 
the result from the fit described in the text. 
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FIG. 4:_Cross section for e+e" -» (a) L>*°L> , (b) _D* + _D" and 
(c) D*D combined. The error bars correspond to statistical 
errors only. 



the D*°D*° channel is estimated by fitting the M 2 ec dis- 
tribution. The fit is performed using a 2™ d -order poly- 
nomial for the background and a signal M r 2 ec lineshape 
obtained from Monte Carlo simulations that reflect the 
composition of the data. The number of ISR candi- 
dates and purities are also summarized in Table III. The 
£)*+£)*- fi na i state has a background consistent with 
zero. 

Due to the small D* + D*~ sample size, the charged 
and neutral mass spectra are summed in Fig. 8. The 
D*D* mass spectrum shows unresolved peaks at -0(4040) 
and 0(4160) and an enhancement at the position of the 
V>(4400) [9]. 

The background is explored using events in the M 2 ec 



sideband regions -2.5 < M r 2 ec < —1.5 GeV 2 /c 4 and 
1.5 < M 2 CC < 2.5 GeV 2 /c 4 , and fitted using Eq. (2). The 
D*D* mass spectrum for these events, normalized from 
the fit to the M 2 CC distribution, is shown as the shaded 
histogram in Fig. 8. 

The D*D* cross section is calculated using the same 
method used to compute the D*D cross section. The 
result, summed over the neutral and charged modes, is 
shown in Fig. 9. All systematic uncertainties which have 
been taken into account for the D*D* mode are listed in 
Table V; the overall uncertainty on the cross section is 
12.4%. 

The D*D* cross section distribution exhibits a thresh- 
old enhancement due to the superposition of the -0(4040) 
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FIG. 5: Am distributions for D*D* candidates after applying 
the |M 2 CC | < 1 GeV 2 /c 4 and m(D*D*) < 6 GeV/c 2 selections, 
for (a) D*° -> D° 7 , (b) D*° -> D°7r°, (c) ZT+ -> D°7r+ with 
D° -» A"-7r+, and (d) ZT+ -» D+tt with D+ -> K~-k + -k + . 
The shaded regions indicate the ranges used to select the D* 
signals. 
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FIG. 7: M r 2 ec distributions for (a) D*°D*° and (b) D* + D*~ . 
The curve in (a) is the result from the fit described in the 
text. 
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FIG. 8: D*D* mass spectrum. The shaded histogram is ob- 
tained from the M 2 sidebands —2.5 < A/ 2 < —1.5 and 



1.5 < Af 2 cc < 2.5 GeV 2 /c 4 . 
fit described in the text. 



The curve is the result from the 



VII. THE DD MASS SPECTRUM 



FIG. 6: M 2 oc distributions for D*°D*° for (a) D*° 



D*° -> D° 7 , (b) D*° 



£°7r°, £° 7 , and (c) £*° 



and V>(4160) resonances. 



In the selection of the D°D° sample we also apply the 
method of resolving ambiguous events having an addi- 
tional ir° and/or 7. Here we veto all events that are 
ambiguous with channels 8-12, obtaining a rejection of 
7.6 % background events in the M r 2 cc | < 1 GeV 2 /c 4 re- 
gion. No such procedure is applied to the D + D~ sample. 
The DD analysis is otherwise identical to that reported 
in Ref. [15]. The resulting Af 2 oc distributions for D°D° 
and D + D~ channels are shown in Fig. 10. The curves 
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TABLE V: Systematic errors, given as fractional errors ex- 
pressed in %, in the evaluation of the D*D* cross section. 



Effect 



Fraction (%) 



Background subtraction 
Branching fractions 
M r 2 c cut 

Particle identification 

Tracking efficiency 

Extraneous tracks 

7T° and 7 reconstruction efficiency 



Total 



2.1 

9.3 
1.3 
2.8 
2.6 
5.7 
1.1 



12.4 
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FIG. 10: M 2 CC distribution for (a) and (b) D+D~ . The 

curves are the results from the fits described in the text. 
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FIG. 9: Cross section e + e" -» D*L>* for combined D*° 1 D*° 
and D* + D*~ . Error bars indicate the statistical errors only. 



are the results from the fits performed using a 2 nd or- 
der polynomial for the background and a M^ ec lineshape 
obtained from Monte Carlo simulations that reflect the 
channel composition of the data. Again, the resulting 
event yields and purities are summarized in Table III. 
The combined DD mass spectrum is shown in Fig. 11. 
The background is explored using events in the A/ r 2 ec side- 
band regions 1.5 < Af 2 cc < 3.5 GeV 2 /c 4 and fitted using 
Eq. (2). This background, normalized from the fit to 
the M r 2 ec distributions, is shown as the shaded histogram 
in Fig. 11. The features in the DD mass spectrum and 
the resulting DD cross section have been extensively dis- 
cussed in our previous publication [15]. 



VIII. FIT TO THE MASS SPECTRA 



Unbinned maximum 
D+D~, D*°D°, D*+D- 



likelihood 



fits to the 
and D*+D* 



D°D°, 
mass 




4 4.5 5 5.5 

m(D D) (GeV/c 2 ) 



FIG. 11: DD mass distribution. The shaded histogram is 
obtained from the M 2 ec sideband. The curve is the result 
from the fit described in the text. 



spectra are performed. We write the likelihood functions 
as 



L = fe 8 {m)\P(m) + ciWi(m)e #1 + c 2V /G'(m)e^ 2 + ... 

+c n W n (m)e^\ 2 + B(m)(l - /), 

(8) 

where m is the DWDM mass, Ci and fa are free param- 
eters, Wi(m) are P-wave relativistic Breit-Wigner distri- 
butions [9], Pirn) represents the nonresonant contribu- 
tion, B(m) describes the background, e B (m) is the av- 
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erage efficiency, and / is the signal fraction fixed to the 
values obtained fitting the M r 2 c distributions. 

The parameters of the 0(4040), -0(4160), and 0(4415) 
are fixed to the values reported in the Review of Particle 
Physics [9]. The parameters of the "0(3770) are fixed to 
the values obtained in our previous analysis of the DD 
system [15]. The DD data require that we include the 
3.9 GeV/c 2 structure, as suggested in Ref. [22], which 
we parametrize empirically as the square root of a Gaus- 
sian times a phase factor J G(m)e 1 ^ 2 . The parameters of 
the Gaussian are fixed to the values obtained in our pre- 
vious analysis of the DD system: m^ggo) = 3943 ± 17 
McV/c 2 , er G(3900) = 52±8 MeV/c 2 [15]. The shape of the 
nonresonant contribution P(m) is unknown; we therefore 
parametrize it in a simple way as 



P(m) = C(m)(a + bm), 



(9) 



where C(m) is the phase space function for £)(*)!)(*), 
and a and b are free parameters. Resolution effects have 
been ignored since the widths of the resonances are much 
larger than the experimental resolution. 

Interference between the resonances and the nonreso- 
nant contribution P(m) is required to obtain a satisfac- 
tory description of the data. The size of the nonresonant 
production is determined by the fit. 

The six D°D a , D+D~ , D*°D°, D*+D-, D*°D*°, and 
D* + D*~ likelihood functions are computed with differ- 
ent thresholds, efficiencies, purities, backgrounds, and 
numbers of contributing resonances appropriate for each 
channel. The fits, summed over the charged and neutral 
final states, are shown in Fig. 12; they provide a good 
description of all the data. In the figure, the shaded ar- 
eas indicate the background estimated by fitting the M r 2 c 
sidebands. The second smooth solid line represents the 
nonresonant contribution where we plot |P(m)| 2 , there- 
fore ignoring the interference effects. The fraction for 
each resonant contribution i is defined by the following 
expression: 



fi 



i\ 2 J \Wi(m)\ 2 dm 



C 3 c k 



J Wj (m) Wj*; (m)dm' 



(10) 



The fractions /, do not necessarily add up to 1 because 
of interference between amplitudes. The error for each 
fraction has been evaluated by propagating the full co- 
variance matrix obtained by the fit. The resulting fit 
fractions and phases are given in Table VI. 



IX. FIT FRACTIONS AND INTEGRATED 
RATES 

The fit gives corrected yields for each charmonium res- 
onance. Since the fits have been performed indepen- 
dently for the neutral and charged modes, the weighted 
means of the fit fractions are used. These can be used 
to compute the integrated rates for each resonance in the 
DD, D*D, and D*D* decay modes, which are reported in 
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FIG. 12: Fits to the (a) DD, (b) D*D, and (c) D*D* mass 
spectra. Data are represented with error bars, the curves 
represent the fitted functions. The shaded histogram corre- 
sponds to the smoothed incoherent background (B(m)) esti- 
mated from sidebands. The second smooth solid line repre- 
sents the nonresonant contribution (|P(m)| 2 ). 



Table VII. The systematic errors take into account un- 
certainties on resonance parameters, Breit-Wigner line- 
shapes, branching fractions, and background estimates. 
The nonresonant contribution has been parametrized in 
an alternative way, P(m) = C(m)e a+bm . Each reso- 
nance parameter has been varied according to its uncer- 
tainty, and the meson radius used in the Blatt-Weisskopf 
damping factor [23], which is present in each relativis- 
tic Breit-Wigner term, has been varied between and 
2.5 GeV -1 . The amounts of backgrounds in the differ- 
ent final states have been varied according to their er- 
rors. The 3.9 GeV/c 2 structure in the DD mass spectrum 
has been alternatively described by a P-wave relativistic 
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TABLE VI: DD, D* D, and D*D* fit fractions (in %) and phases. Errors are statistical only. 



DD DD D* D D*D D*D* D*D* 

Res. fraction phase fraction phase fraction phase 

\P\ 2 38.5 ± 7.1 0. 49.9 ± 5.6 0. 56.8 ± 9.2 0. 

V>(3770) 31.3 ± 3.3 1.58 ± 0.46 

G(3900) 23.9 ± 5.8 5.46 ± 0.64 

0(4040) 31.2 ± 5.3 1.39 ± 0.55 34.5 ± 6.0 1.74 ± 0.33 5.7 ± 4.4 3.37 ± 0.48 

0(4160) 3.1 ± 3.3 2.75 ± 0.58 12.2 ± 3.8 2.26 ± 0.16 30.6 ± 7.3 5.94 ± 0.33 

0(4400) 2.0 ± 1.2 3.38 ± 0.37 0.6 ± 0.7 4.37 ± 0.47 3.6 ± 2.4 5.03 ± 0.45 

Sum 130 ± 12 97 ± 8 97 ± 13 



TABLE VII: Integrated rates (in nb- MeV) for e+e~ -> 0(4040), e+e~ 0(4160), and e+e" 0(4400) in the DD, D*D, and 
D*D* decay modes. The first error is statistical, the second systematic. 



Decay mode 


0(4040) 


0(4160) 


0(4400) 


DD 


11.0 ± 1.8 ± 5.6 


1.0 ± 1.3 ± 1.0 


0.5 ±0.3 ±0.1 


D*D 


46.6 ±7.0 ±4.9 


13.8 ± 4.4 ± 1.5 


0.6 ±0.8 ±0.1 


D*D* 


8.3 ±6.4 ±1.0 


40.6 ± 9.7 ± 5.0 


3.6 ±2.4 ±0.4 



Breit Wigner with free parameters. This effect domi- 
nates the systematic uncertainty on the 0(4040) rate in 
the DD mass spectrum. The deviations from the central 
value are added in quadrature. Systematic effects also 
include the uncertainty on the total cross sections. 

The corrected yields can also be used to compute the 
branching fraction ratios. The results are shown in Ta- 
ble VIII together with predictions of models: signifi- 
cant discrepancies are observed, expecially with the 3 Po 
model [5]. 



ment from the J/ipn + TT channel [2]. We obtain 



LIMITS ON THE DEC AYS_Y (4260) 
AND U(4260) -> D*D* 



D*D 



The D*D and D*D* mass spectra have been refit with 
an additional Y(4260) resonance, which is allowed to in- 
terfere with all the other terms. 

The fit gives F(4260) fractions of (2.2 ± 2.9 sta t ± 
2.5 syst )% and (4.0 ± 2.0 sta t ± 4.2 syst )% corresponding 
to 18 ± 24 sta t ± 21 syst and 9 ± 5 sta t ± 10 sys t events for 
Y(4260) -> D*D and Y(4260) -> D*D* , respectively. 
Systematic errors due to uncertainties on masses and 
widths of the 0(4040), 0(4160), 0(4415), and F(4260) 
resonances are evaluated by varying the masses and 
widths by their uncertainty in the fit. The amount of 
background in each final state is varied within its statis- 
tical error, and the meson radii in Breit- Wigner terms 
are varied between and 2.5 GeV" 1 . Deviations from 
the central value are added in quadrature. 

These F(4260) yields in the D*D and D*D* channels 
are used to compute the cross section times branching 
fraction, which can then be compared to our measure- 



and 



ff(y(4260) -> D*D) 
B(Y (4260) -> J/ipir+iT-) 



£(Y(4260) -» D*D*) 
S(F(4260) -» J/ipn+TT-) 



< 34, 



< 40, 



(11) 



(12) 



at the 90% confidence level. 

Using the DD cross section measured in the earlier 
BABAR work [15], we obtain the sum of the e + e~ — > DD, 
e + e~ — > D*D, ande + e~ — > D*D* cross sections shown in 
Fig. 13: the arrow indicates the position of the F(4260), 
which falls in a local minimum, in agreement with the 
cross section measured for hadron production in e + e _ 
annihilation [9]. 



XI. CONCLUSIONS 

We have studied the exclusive ISR production of the 
DD, D*D, and D*D* systems. The mass spectra show 
production of the J PC = 1~ states -0(3770), V(4040), 
■0(4160), and -0(4415). Fits to the mass spectra pro- 
vide amplitudes and relative phases for the charmo- 
nium states, from which first measurements of branch- 
ing fraction ratios are obtained. Finally, upper limits 
on F(4260) -> D*D and F(4260) -> D*D* decays are 
computed. 

If the y(4260) is a 1 charmonium state, it should de- 
cay predominantly to DD, D*D, and D*D* [5, 6]. Within 
the present limited data sample size, no evidence is found 
for F(4260) decays to DD, D*D, or D*D* . Other expla- 
nations for the F(4260) have been proposed, such as a 
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TABLE VIII: Ratios of branching fractions for the three tfi resonances. The first error is statistical, the second systematic. 
Theoretical expectations are from the 3 Po model [5], C 3 model [6], and pKp model [14]. 

Ratio measurement 3 Po C 3 and pKp 

1) £(V(4040) DD)_/B{ip (4040) -»■ D*D)_ 0.24 ± 0.05 ± 0.12 0.003 0.14 [14] 

2) S(V(4040) -> L»*D*)/B(V'(4040) -> ITD) 0.18 ± 0.14 ± 0.03 1.0 0.29 [14] 

3) S(V>(4160) -> DD)_[B(ip(4\Q0) ->■ D*D^) 0.02 ± 0.03 ± 0.02 0.46 0.08 [6] 

4) Z3(V>(4160) -> D^D)/B(i/>(4160) ->■ D*D*) 0.34 ± 0.14 ± 0.05 0.011 0.16 [6] 

5) B(V>(4400) — ► DD)JB{i '(4400) -> D*^_) 0.14 ± 0.12 ± 0.03 0.025 

6) g(V>(4400) -> P*P)/g(V>(4400) -» D*D*) 0.17 ± 0.25 ± 0.03 0.14 



hybrid, baryonium, molecule or tetraquark state. In the 
case of a hybrid state, the decay rates to DD, D*D, and 
D*D* are expected to be small [10, 24]. 
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FIG. 13: Sum of e+e" -> DD, e + e~ -> L>*D, and e + e~ -> 
D*D* cross sections. The arrow indicates the position of the 
F(4260). 



XII. ACKNOWLEDGEMENTS 

We are grateful for the extraordinary contributions of 
our PEP-II colleagues in achieving the excellent luminos- 
ity and machine conditions that have made this work pos- 
sible. The success of this project also relies critically on 
the expertise and dedication of the computing organiza- 
tions that support BaBAR. The collaborating institutions 
wish to thank SLAC for its support and the kind hospi- 
tality extended to them. This work is supported by the 
US Department of Energy and National Science Foun- 
dation, the Natural Sciences and Engineering Research 
Council (Canada), the Commissariat a PEnergie Atom- 
ique and Institut National de Physique Nucleaire et de 
Physique des Particules (France), the Bundesministerium 
fiir Bildung und Forschung and Deutsche Forschungsgc- 
mcinschaft (Germany), the Istituto Nazionalc di Fisica 
Nucleare (Italy), the Foundation for Fundamental Re- 
search on Matter (The Netherlands), the Research Coun- 
cil of Norway, the Ministry of Education and Science of 
the Russian Federation, Ministerio de Education y Cien- 
cia (Spain), and the Science and Technology Facilities 
Council (United Kingdom). Individuals have received 
support from the Marie-Curie IEF program (European 
Union) and the A. P. Sloan Foundation. 



[1] S.-K. Choi et ai, Belle Collaboration, Phys. Rev. Lett. 

91, 262001 (2003). 
[2] B. Aubert et ai, BABAR Collaboration, Phys. Rev. Lett. 

95, 142001 (2005). 
[3] B. Aubert et ai, BABAR Collaboration, Phys. Rev. Lett. 

98, 212001 (2007). 

[4] X.L. Wang et al., Belle Collaboration, Phys. Rev. Lett. 

99, 142002 (2007). 

[5] T. Barnes, S. Godfrey, and E.S. Swanson, Phys. Rev. 

D72, 054026 (2005). 
[6] E.J. Eichten, K. Lane, and C. Quigg, Phys. Rev. D73, 



014014 (2006). 

[7] Charge conjugate states are implied throughout this 
work. 

[8] X.H. Mo et al, Phys. Lett. B640, 182 (2006). 
[9] C. Amsler et al., Review of Particle Physics, Phys. Lett. 
B667, 1 (2008). 
[10] S. L. Zhu, Phys. Lett. B625, 212 (2005); E. Kou and 
O. Pene, Phys. Lett. B631, 164 (2005); F. E. Close and 
P. R. Page, Phys. Lett. B628, 215 (2005). 
[11] C.F. Qiao, Phys. Lett. B639, 263 (2006). 
[12] L. Maiani, F. Piccinini, A. D. Polosa, and V. Riquer, 



16 



Phys. Rev. D72, 031502 (2005). 
[13] X. Liu, X. Q. Zeng, and X. Q. Li, Phys. Rev. D72, 054023 
(2005). 

[14] E.S. Swanson, Physics Reports 429, 243 (2006). 
[15] B. Aubert et al, BABAR Collaboration, Phys. Rev. D76, 
111105 (2007). 

[16] G. Pakhlova et al, Belle Collaboration, Phys. Rev. D77, 
011103 (2008). 

[17] G. Pakhlova et al, Belle Collaboration, Phys. Rev. Lett. 

98, 092001 (2007). 
[18] J. Libby et al, CLEO Collaboration, Nucl. Phys. Proc. 



Suppl. 181-182, 127 (2008). 

[19] B. Aubert et al, BABAR Collaboration, Nucl. Instrum. 
Methods A479, 1 (2002). 

[20] S. Agostinelli et al, GEANT Collaboration, Nucl. In- 
strum. Methods Phys. Res., Sect. A 506, 250 (2003). 

[21] M. Benayoun et al, Mod. Phys. Lett. A14, 2605 (1999). 

[22] E. Eichten et al, Phys. Rev. D21, 203 (1980). 

[23] J.M. Blatt and W.F. Weisskopf, Theoretical Nuclear 
Physics, John Wiley & Sons, New York, 1952. 

[24] F. Iddir, L. Semlala, hep-ph/0611183. 



